Epitaxial Growth Kinetics with Interacting Coherent Islands 
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The Stranski-Krastanov growth kinetics of undislocated (coherent) 3-dimensionaI islands is studied 
with a self-consistent mean field rate theory that takes account of elastic interactions between the 
islands. The latter are presumed to facilitate the detachment of atoms from the islands with a 
consequent decrease in their average size. Semi-quantitative agreement with experiment is found 
for the time evolution of the total island density and the mean island size. When combined with 
scaling ideas, these results provide a natural way to understand the often-observed initial increase 
and subsequent decrease in the width of the coherent island size distribution. 
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Heteroepitaxy begins with the formation of a thin, 
lattice-matched wetting layer if the energy gain from 
substrate-adlayer adhesion exceeds the elastic energy cost 
from lattice constant misfit Sa/a. As deposition contin- 
ues, 2-dimensional (2D) islands nucleate on top of the 
wetting layer. These islands contribute to the build-up 
of elastic strain and, for this reason, the system does 
not tolerate their growth, coalescence, and re-nucleation 
indefinitely. Instead, at large misfit, coherent (undislo- 
cated) 3-dimensional (31?) islands form that are lattice 
matched near their base but are largely strain-relieved 
near their top and sidewalls. Further deposition leads 
to their growth and eventual coalescencei- This is the 
so-called Stranski-Krastanov growth modetl. 

A coherent island is the source of strain fields because 
it elastically distorts the wetting layer and substrate in 
its immediate vicinity. Early on, several experimental 
groups observed a significant decrease in the mean size 
of the coherent islands at relatively early stages of growth 
and suggested the possible role of long, range strain fields. 
For example, Ponchet and co-workersB presented data for 
the InAs/InP(001) system and pointed out that elastic 
interactions should cause islands to destabilize one an- 
other because their interactions are mutually repulsive. 
Kobayashi and co-workersu identified several other fea- 
tures of the island-island interaction as a basis for un- 
derstanding their experiments on the InAs/GaAs(001) 
system. Theoretical work on island interactions has been 
limited to equilibrium considerationsu up to the present 
time. 

In this paper, we present a theoretical analysis of 
Stranski-Krastanov growth kinetics that generalizes pre- 
vious work by Dobbs and co-workersO to take acccount 
of island interactions and atom detachment from 31? is- 
lands. Dobbs et al. employed a mean field theory for 
the density of adatoms ria, the density of 2D islands n2, 
their average size S2, the density of SD coherent islands 
n3, and their average size S3. A rate equation was derived 
for each based on the physical processes of adatom depo- 
sition, surface diffusion, attachment and detachment of 



adatoms from the islands, etc. In brief, an incident flux 
F contributes directly to the increase of adatom popua- 
tion. The adatoms diffuse on the surface with a diffusion 
constant D = ojo? e'xjp{—Es/kBT) where uj is an attempt 
frequency, a is the lattice constant, Es is the energy bar- 
rier for diffusion, ks is the Boltzmann constant and T is 
the temperature. Diffusing adatoms that meet bond to- 
gether to form small 2D islands but thermal fluctuations 
can cause them to break apart if the island size is too 
small. There is a critical island size i such that islands 
of size i and less are unstable. 

An island grows by capturing adatoms from both the 
vapor and the substrate. The rates for these processes are 
Fk and Dana , where k is the direct capture number and 
a is the diffusion capture number. To relieve strain, 2D 
islands convert into 31? islands at a rate 72. We assume 
that atoms that detach from the edges of a 2D island 
do not leave the island but instead migrate to the top of 
the island. On the other hand, we suppose that atoms do 
detach from 31? islands (at a rate I/T3) when interactions 
become significant. A fraction 7712 of these attach to 2Z? 
islands. The remaining fraction nia contributes to the 
adatom population. In this work, we approximate 7712 by 
the areal coverage of 21? islands. 

Rate equations that incorporate all of these elementary 
processes are 
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The suffixes i, 2 and 3 for n and a denotes critical nuclei 
(of size i), 2D and 31? islands respectively. 

We assume that 2D islands are circular with radius r 
and 3D islands are truncated pyramids with base length 
Z, height h an d base angle (j). The radius of a 2D island 
is r = V'sT/Tr. We assume that ZD islands very quickly 
achieve their equilibruim shape and that the angle (j) does 
not change significantly during growth. For a given is- 
land size S3, h and I are found by minimizing the energy 
expression jderived for a ZD coherent island by Tersoff 
and TrompQ. 

The direct capture number k is given by the surface 
area of the island normal to the incident flux, i.e., K2 = S2 
and K3 = l'^. The diffusion capture number a measures 
the efficiency with which an island captures adatoms 
from the surface. We coHipute a using the prescription 
of Bales and co-workerst2l which relates it to the size of 
the diffusional depletion zone ^ that surrounds each is- 
land. For a circular island of radius r with no barrier to 
adatom attachment we solve self-consistently 



cr = 27r 
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where K^ix) is the modified bessel function of order n. 
We use equation (g) for circular 2D islands and for 173 as 
well (with r replaced by Z/2) because the details of the 
island shape should not affect the results significantly. 

Conversion of a 2D island to a 3D island occurs when 
a sufficient density of atoms is present on its top (due to 
strain-driven detachment from its perimeter and upward 
migration) to nucleate a lapw island at its center. The 
requisite nuclcation rate id 



72 = ^r2Dexp[(D, - (i + \)Ei{r))l(kBT)\ 
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where Ei is the binding energy of critical nuclei and 
E^(r) ~ Eq \u{r / a) I (r / a) is a size dependent energy 
barrier for the detachment of atoms from the 2D island. 
The form of 72 as a function of r is such that 2D islands 
barely convert at all until they reach a size s* after which 
most of them convert very rapidly. 

The escape rate of an atom from a 3D island is 



i = :^ex ( ^^ ) 



(5) 



where Ef, is .-the energy barrier for detachment. Elas- 
ticity theoryO predicts that the change in the barrier 
due to strain is AD;, w (ry„ — ?7s)e, where e is the local 
strain and r]s{T]u) is the local surface stress at the binding 
site (transition state) configuration. The predicted lin- 
earity with jStrain has been confirmed by first principles 
calculationsu. The strain field due to a misfitting island 
is proportional to the size of the island and varies as d~^ 
for distances d far from the islandQ. We therefore put 



where E^ is the strain independent part of the energy 
barrier and ds = l/y^7i3 is the average 3D island sepa- 
ration. We treat a as an adjustable parameter because 
the surface stress difference discussed above is difficult to 
estimate. The factor of 27r is, in this model, the mean 
number of islands that are nearest-neighbors to a given 
island. 



10 



10" 



10 



10° 



10' 



10 



10 



(a) 



— 1 1 — 

.05 .10 

e(ML) 



— 1 — 
.15 



10 



(b) 



1.55 



1.60 



1.65 
e(ML) 



1.70 



1.75 



FIG. 1. 3D Island Densities:(a) from the present theory, 
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(b) from the data of Kobayashi et. al 



The rate equations were integrated numerically us- 
ing an aljgorithm suited for systems of stiff differential 
equationalj. We used values of the parameters typical of 
those found in experiments, T = 900 K, i^ = 0.1 ML/s, 
a = 3.0 A, i = 4, Do = 3.5 ev, E, = 0.5 ev, D, = 1.0 
ev, DO = 0.7 ev, </> = 25° and Sa/a = 0.05. Our resuhs 
for the time (coverage) evolution of the 3D island density 
and mean size are shown in Figure ^(a) and Figure 0(a). 
For comparision we have plotted the experimental results 
of Kobayashi et al. for InAs/GaAs(001)l in Figure |l|(b) 
and 0(b). The sizes were estimated from the published 
experimental distributions of island heights and island 
widths. Note also that we have shifted the theoretical 
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FIG. 2. Average 3D Island Size: (a) from the present the- 
ory, (b) estimated from the data of Kobayashi et. aB 



curves to align the rapid island density onsets, because 
the precise onset position is related to alloying)!^ that we 
do not attempt to model. 

The 3D island density initially rises very rapidly due 
to the fast conversion of 2D islands to 31? islands. It 
then tends to saturate because, as a result of conversion, 
the average 2D island size decreases below s* . During 
this time the average 3D island size continue to grow. 
Soon the interactions become important and significant 
detachment of atoms from the 3D islands begins. This 
results in the very rapid decrease of S3 seen in Figure 
H(a). The detached adatoms that re-attach to 2D islands 
increase the average size of the latter to s* which, in 
turn, leads to more 2D to 31? conversion. That is why 
the 3D island density increases again. The same trend is 
seen in the experimental data although we do not obtain 
quantitative agreement between our model and the data. 

The results shown correspond to a ~ 120 ev which 
is three orders of magnitude greater than typical elastic 
energies. This large number arises in our model because 
the rapid decrease, in 3D island size seen in the data 
of Kobayashi et alB occurs when the experimental mean 
island separation is ten times larger than the mean island 



radius! Of course, the real system has many islands at 
much closer distances than our simple mean theory can 
describe, but it remains the case that detachment effects 
seem to set in far earlier than simple elasticity estimates 
would suggest. The detailed origin of this behavior is an 
outstanding open question and our simple form (|6|) must 
be regarded as a convenient parameterization. 

In principle, the entire island size distribution can be 
gotten from a rate equation analysis. In practice however, 
it is prohibitively difficult to solve the tens of thousands 
of equations so generated. This theoretical problem is 
ameliorated for the case of 2D homoepitaxy because the 
island size distribution shows scaling behaviorli3. It is 
therefore highly signficant that Ebiko et. al have shown 
that the 3D coherent island size distribution for the 
InAs/GaAs(001) system also shows scalingEJ. Their data 
fits remarkably well to,an analytic scaling form suggested 
for 2D homoepitaxytJ. In detail, the number of islands 
of size s, n, takes the form 



(s)^ (s) 
where Oc — '^^ sns and 

f{u) = l.luexp(-0.27u3-^). 
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It is surprising that an island distribution that works well 
for 2D islands works equally well for 3D islands. Even 
more puzzling is the fact that (@) applies only to situa- 
tions where atom detachment from 2D islands is strictly 
forbidden {i — 1) whereas the coherent islands studied 
here shrink precisely due to copious detachment. 

This can be understood if we parameterize the island 
size distribution not by a fictitious "critical island size" 
but by the ratio of the net detachment rate from an island 
to the net attachment rate to an islandlij, namely 



A = 



1/7 



(Fk + Daria)' 
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Monte Carlo simulations of 2D homomepitaxy show 
that A jiarameterizes a continuous family of scaling 
functionsta. When A ~ 1 or less, the island size dis- 
tribution fits ([7|) and (JSJ) very well even when significant 
detachment is present. The computed time evolution of 
A for our model is shown in Figure y. Note that its value 
exceeds unity when the island interactions are most im- 
portant but only barely so. This is not inconsistent with 
the rapid decrease in the average coherent island size seen 
in Figure (|2|) because the rate equation for this quan- 
tity in (1) involves the difference (rather than the ratio) 
of the attachment and detachment rates which can be 
large. These considerations provide jSu rationale for the 
fitting procedure used by Ebiko et a/.ll3. 

We conclude that island interactions strongly influence 
the average island size but not the island size distribution 
scaling function. This is important because it means that 
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FIG. 3. Coverage dependence of A, the detachment rate 
to attachment rate ratio for a 3D island. 
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we can " synthesize" the time dependence of the entire is- 
land size distribution merely from knowledge of the time 
dependence of the average size. This is shown in Figure 
(Q). As expected, the island size distribution broadens 
and its peak position moves to the right as the coverage 
increases from zero. But as a consequence of equations 
(^ and (P) , the decrease of S3 when interactions become 
imporant induces a narrowing of the distribution and a 
shift back to the left. Precisely this behaiaor is seen in 
the experimental island size distributionaa'cl. 

In summary, we have generalized the theory of Dobbs 
et. aln to take account of island-island elastic interac- 
tions that are presumed to induce atom detachment from 
3D coherent islands. Semi-quantitative agreement was 
found with experimental results for InAs/GaAs(001) but 
the large value for the interaction parameter needed to 
model the data suggests that we still lack a good under- 
standing of the energy barriers to detachment for this 
problem. In conjunction with a scaling ansatz, the re- 
sults could be used nonetheless to rationalize the ubiqui- 
tous " narrowing" of the full island size distribution seen 
in experiment. An interesting and open question is to 
establish the veracity of this scaling assumption in a the- 
oretical framework. 
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